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(Dated: November 4, 2008) 

We report results of a search for the pair production of the hghtest supersymmetric partner of 
the top quark, ii, using a data set corresponding to an integrated luminosity of 1 flj^^ collected 
by the DO detector at a pp center-of-mass energy of 1.96 TeV at the Ferniilab Tevatron Collider. 
Both scalar top quarks are assumed to decay into a 6 quark, a charged lepton and a scalar neutrino. 
The search is performed in the electron plus muon and dielectron final states. The signal topology 
consists of two isolated leptons, missing transverse energy, and jets. We find no evidence for this 
process and exclude regions of parameter space in the framework of the minimal supersymmetric 
standard model. 
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PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 



Supersymmetric theories [ij predict for every standard 
model (SM) particle the existence of a superpartncr that 
differs by half a unit of spin. The top quark would have 
two scalar partners, and i^, corresponding to its left- 
and right-handed states. Mixing between II and in, be- 
ing proportional to the top quark mass mt, may lead to a 
possible large mass splitting between the physical states 
ti and Hence, the lightest supersymmetric partner of 
the top quark, ti, might be light enough to be produced 
at the Fermilab Tevatron collider. 

In this Letter we present a search for scalar top (stop) 
pair production in a data sample corresponding to an 
integrated luminosity of 1 fb~^ collected at a centcr-of- 
mass energy of 1.96 TeV with the DO detector during 
Run II of the Fermilab Tevatron pp collider. The phe- 
nomenological framework is the minimal supersymmet- 
ric standard model (MSSM) with R-parity conservation. 
We assume that BR{ti biv) = 1, where v is 

the scalar neutrino (sneutrino). Among possible stop de- 
cays Q, this final state is one of the most attractive; in 
addition to a 5 quark, it benefits from the presence of a 
lepton with high transverse momentum with respect to 
the beam axis {pr)- The sneutrino is either the light- 
est supersymmetric particle (LSP) or decays invisibly: 
V — > vx\ or vG where the lightest neutralino, x^, or the 
gravitino, G, is the LSP. We suppose an equal sharing 
among lepton flavors and consider titi — > bhU vv fi- 
nal states, with II = e^/z^ (e/i channel) and II = e+e" 
(ee channel) . The signal topology consists of two isolated 
leptons, missing transverse energy (I^t), coming mainly 
from undetected sneutrinos, and jets. A search for stop 
pair production in the and /i/.t (tiii — > bbfifii'i') chan- 
nels has previously been performed by the DO collabora- 
tion using a data set corresponding to a luminosity of 
428 pb""'^. The e/i sample in |3i| is a subset of the data 
sample used in this analysis. Searches for stop pair pro- 
duction in the bb££ vv final state have been reported by 
the ALEPH, L3, and OPAL coUaborations 

The DO detector [sj comprises a central tracking sys- 
tem surrounded by a liquid- argon/uranium sampling 
calorimeter and muon detectors. Charged particles 
are reconstructed using multi-layer silicon detectors and 
eight double layers of scintillating fibers in a 2 T mag- 
netic field produced by a superconducting solenoid. Af- 
ter passing through the calorimeter, muons are detected 
in the muon system comprising three layers of tracking 
detectors and scintillation counters. Events containing 
electrons or muons are selected for offline analysis by an 
online trigger system. A combination of single electron 
(ee channel) and dilepton (e^ channel) triggers is used to 
tag the presence of electrons and muons based on their 
energy deposition in the calorimeter, hits in the muon 
detectors, and tracks in the tracking system. 



In pp collisions, stops are pair-produced via quark- 
antiquark annihilation and gluon fusion. The ii pair pro- 
duction cross section, cr^^j^, depends primarily on m^^, 
with only a weak dependence on other MSSM parame- 
ters. At \/s = 1.96 TeV, err ~ at next-to-leading-order 
(NLO), calculated with prospino [6|, ranges from 15 pb 
to 0.5 pb for 100 < m^^ < 180 GeV. These cross sec- 
tions are estimated using CTEQ6.1M parton distribution 
functions (PDF) 0,11] and equal renormalization and fac- 
torization scales p,rj = rni_^ . A theoretical uncertainty of 
about 18% is estimated due to scale and PDF choice. 

Three-body decays of the stop are simulated using 
COMPHEP ^ and PYTHIA [l^] for parton-level generation 
and hadronization, respectively. We consider a range of 
stop mass values from 100 to 200 GeV in steps of 10 GeV. 
The range of sneutrino masses explored extends from 40 
to 140 GeV in steps of 10 to 20 GeV. For each choice of 
, 10,000 events are generated. Background pro- 



'ti' 



cesses are simulated using the pythia and ALPGEN [11 | 
Monte Carlo (MC) generators. ALPGEN is interfaced with 
PYTHIA for parton showering and hadronization. The 
MC samples use the CTEQ6L PDF and are normalized 
using next-to- leading order cross sections 13, 13, IM ■ 
generated events are passed through the full simulation of 
the detector geometry and response based on GEANT fl^ . 
MC events are then reconstructed and analyzed with the 
same software as used for the data. 

The signal topology depends both on and on 
the mass difference Am = mi_^ — mt,- The pT of the 
leptons and b quarks decrease with smaller values of 
Am and I^t values are correlated with mi_^ and Am. 
For both e/i and ee channels, the two signal points 
[mi^, mi,] = (140,110) GeV and (170,90) GeV, referred 
to respectively as "Signal A" and "Signal B" in the fol- 
lowing, are chosen to illustrate the effect of the selections 
for low TO(^ and low Am (Signal A) and for high m^^ and 
high Am (Signal B). 

The main SM background processes mimicking the sig- 
nal signature are Z/^* t^t^, WW, WZ, ZZ, and 
tt (e/i and ee decay channels), Z/j* — > e~'"e~ (ee chan- 
nel), and instrumental background (e/i and ee channels). 
All but the latter are estimated using MC simulations. 

Electrons are identified as clusters of en- 
ergy in calorimeter cells in a cone of size 
n = v/(A(?!))2 + (A77)2 = 0.4 where^is the az- 
imuthal angle and r/ the pseudorapidity [l6|. Electron 
candidates are required to have a large fraction of their 
energy deposited in the electromagnetic layers of the 
calorimeter. The clusters are required to be isolated 
from hadronic energy depositions. The calorimeter 
isolation variable / = [£'tot(0.4) - £;em(0.2)]/£:em(0.2) 
is less than 0.15, where _Etot(0.4) is the total transverse 
shower energy in a cone of radius TZ = 0.4 and £'em(0.2) 
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FIG. 1: Distributions {ejj, channel) of (a) I^t after preselection, (b) A(f){e,l^T) and (c) A(^(^, _^t) after Emu 1, (d) 
A4i{e, I^T)+A(j){^, I^t) after Emu 2, (e) Ht and (f) St after Emu 3, for observed events (dots), expected background (filled 
areas), and signal expectations for Signal A (solid line) and Signal B (dashed line). 



is the electromagnetic energy in a cone TZ = 0.2. The 
clusters are also required to have a spatially-matched 
track in the central tracking system with pT larger than 
8 GeV, and to have a shower shape consistent with that 
of an electron. Electrons are also required to satisfy 
identification criteria combined in a likelihood variable 
and based on multivariate discriminators derived from 
calorimeter shower shape and track variables. Only 
central electrons {\r]\ < 1.1) with transverse energy 
with respect to the beam axis {Et) measured in the 
calorimeter larger than 15 GeV are considered. 

Muons are reconstructed by finding tracks pointing to 
hit patterns in the muon system. Non-isolated muons are 
rejected by requiring the sum of the transverse momenta 
of tracks inside a cone of radius 7^ = 0.5 around the 
muon direction to be less than 4 GeV, and the sum of 
transverse energy in the calorimeter in a hollow cone of 
size 0.1 < TZ < 0.4 around the muon to be less than 
4 GeV. To reject cosmic ray muons, requirements on the 
time of arrival of the muon at the various scintillator 



layers in the muon system are made. Muons with |?7| < 2 
and pt > 8 GeV are considered. 

Jets are reconstructed from the energy de- 
position in the calorimeter towers using the 
Run II cone algorithm [l7| with a radius 
7^cono = x/(A0)2 + (Ay)2 = 0.5, where y is the 
rapidity 16]. Jet energies are calibrated to the particle 
level using correction factors primarily derived from 
the transverse momentum balance in photon plus jets 
events. Only jets with pT > 15 GeV and \ri\ < 2.5 are 
considered. The I^t is calculated using all calorimeter 
cells and is corrected for the jet and electromagnetic 
energy scales and for the momentum of selected muons. 

In each event, the best primary vertex is selected 
from all reconstructed primary vertices as the one with 
the smallest probability of originating from a mini- 
mum bias interaction [18]. Its longitudinal position 
with respect to the detector center, z, is restricted to 
\z\ < 60 cm to ensure efficient reconstruction. The lep- 
tons in an event are required to be isolated from each 
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other (7e(f,f) > 0.5) and from a jet (7e(f, jet) > 0.5). 

The instrumental background is due to either misiden- 
tified electrons or muons, mismeasured I^Ti or electrons 
or muons from multijet processes that pass the lepton 
isolation requirements presented above. Data samples 
dominated by instrumental background are selected by 
inverting the muon isolation requirements or the electron- 
likelihood cut (e/j, channel) or both electron-likelihood 
criteria (ee channel). The normalization factors for those 
samples are estimated from observed events. In the 
e/Lt channel, an exponential fit is performed to the I^t dis- 
tribution in the range If^x < 35 GeV, after subtraction of 
the MC estimates of the non-instrumental backgrounds, 
in events containing one electron and one muon. In the 
ee channel, the normalization is performed using both 
electron Et shapes in events containing two electrons in 
a domain where the instrumental background has a large 
contribution. 

The integrated luminosity [3] of the e/j, data sample is 
1100 ± 67 pb~^. Events are preselected with the require- 
ment that they contain one electron and one muon. To 
remove a large part of the instrumental background as 
well as events coming from Z/^* t^t~ , selections on 
the If^T [Fig- [ija)] and on the I^t significance, S{I^t), 
defined as the ratio of the I^t in an event to its esti- 
mated uncertainty given the expected resolutions on the 
Pt measurements for the selected leptons and jets, are 
applied: 

> 30 GeV 
S{^t) > 4. (Emu 1) 

At this stage, the instrumental and 
Z/Y T^T^ events comprise a large part (41%) 
of the total background. In these processes, recon- 
structed leptons are correlated with the I^t, giving rise 
to higher event populations at high and low values of 
the azimuthal angle difference between the leptons and 
I^T, with a low value of the angular difference for one 
lepton being correlated with a high value for the other. 
As there is a higher background contribution at low 
values of the angular distributions [Figs. [T]^b) and[TJc)], 
we require: 

A(j){ii, ^t) > 0.4 rad 

A(j){e, ^t) > 0.4 rad. (Emu 2) 

To reduce the Z/j* t^t^ background, selections on 
the transverse mass of the muon and ^t, Mt{ij., Ifr) 0|j 
and of the electron and Ifx: Mx^e, I^t), are applied. To 
further reduce this background, we use the azimuthal an- 
gular differences between the leptons and the missing en- 
ergy, A0(/i, ^t) and A(j){e, I^t), which should be large 
[Fig. md)]. We require: 



Mt{h,$t) > 20 GeV 
MT{e,^T) > 20 GeV (Emu 3) 
A(/)(/i,^T) + A0(e,^T) > 2.9 rad. 

The number of events surviving at each analysis step 
for the data, for each background component, and for 
the two signal samples A and B are summarized in Ta- 
ble |T1 After all selections, the WW, tt, and instrumen- 
tal background contributions dominate. To separate the 
signal from these backgrounds, two topological variables 
are used: St, defined as the scalar sum of the muon 
Pt, the electron pt, and the I^t', and Ht, defined as 
the scalar sum of the transverse momenta of all the jets. 
WW and instrumental backgrounds populate low values 
of Ht and St while top quark pairs have large values 
for both variables. The signal distribution depends on 
the stop mass and on the mass difference Am, with low 
values of Am having low values of Ht and St [Figs. [Ije) 
and [Iff)]. Rather than selecting events using these two 
variables, the numbers of events predicted for signal and 
background are compared to the observed numbers in 
twelve [St,Ht] bins (Table HI)) when extracting limits on 
the cross section for the e/i channel. 

The integrated luminosity of the ee data sample is 
1043 ± 64 pb~^. At preselection, two electrons are re- 
quired. Z/'j* e+e~ events account for 94% of the 
total background. While the signal is characterized by 
the presence of jets originating from the hadronization of 
b quarks, the Z/^* e^e~ background owes the pres- 
ence of jets to gluons from initial state radiation which 
hadronize into softer jets, resulting in a lower multiplic- 
ity of jets. To keep sensivity to low Am signals while 
rejecting substantial background, we require at least one 
jet [Fig. Ha)]: 

7V(jets) > 1. (Dielec 1) 

To reject contributions from both the instrumental and 
Z/7* — ^ e+e^ backgrounds, cuts on the I^t and on its 
significance are performed: 

> 15 GeV 
S{^t) > 5. (Dielec 2) 

At this stage of the analysis, the Z/j* e+e" sample is 
still dominant [Fig. Hfb)] and give rise to higher event 
populations at high values of the azimuthal angle differ- 
ence between the two electrons. To remove these events, 
the following selection is applied: 

A(j){ee) < 3 rad. (Dielec 3) 

To increase the search sensitivity in this channel, we 
take advantage of the presence of jets originating from 
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TABLE I: Numbers of events observed in data and expected from SM background processes and the two signal samples A and 
B at the various stages of the analysis in the e/i channel. The quoted uncertainties are statistical only. 







Total SM 


Backf 


ground contributions 








Selection 


Data Background ^/7* 


— * t'^t" 


tt Diboson Instrumental 


Signal A 


Signal B 


Preselection 


735 


736 ± 15 


458 


29.7 60.6 


188 


34.0 ± 1 


26.3 ±0.7 


Emu 1 


106 


106 ±5 


23 


23.5 38.7 


21 


10.6 ±0.7 


19.4 ±0.6 


Emu 2 


71 


77 ±4 


5.9 


20.0 36.2 


15 


8.4 ±0.7 


17.6 ±0.6 


Emu 3 


61 


65 ±4 


0.7 


16.4 34.5 


13 


6.0 ±0.6 


16.1±0.5 
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FIG. 2: Distributions (ee channel) of (a) the number of jets after the preselection, (b) Ifx after Dielec 2, (c) the dielectron 
invariant mass and (d) Ifr after Dielec 3, (e) Ht and (f) St after Dielec 5, for observed events (dots), expected background 
(filled areas), and signal expectations for Signal A (solid line) and Signal B (dashed line). 



the fragmentation of long-lived h quarks in the signal. A 
neural network (NN) tagging tool [21|] for heavy flavor 
that combines information from several lifetime-based fa- 
taggers to maximize the h quark tagging efRciency is used 
for this purpose. At least one jet in the event is required 
to be 6-tagged (Dielec 4) by satisfying a given NN selec- 
tion. The b quark tagging operating point preserves high 
efficiency for the detection of b jets (« 66%) with a « 3% 
probability for a light parton jet to be mistakenly tagged. 
This point maximizes the sensitivity of the analysis for 
stop masses of 130 to 140 GeV and for low Am. At this 



stage, most of the surviving Z/^* e+e~ events have a 
dielectron mass in the vicinity of the Z boson resonance 
and low I^t values [Figs, dfc) and^d)]. To further sup- 
press this background while preserving the signal, a cut 
in the plane [M{e,e),l^T] is applied. This selection is 
optimized for low Am signals and is defined by: 

M(e, e) ^ [75, 105] GeV if < 30 GeV. (Dielec 5) 

The selections applied in the ee channel are summa- 
rized in Table IIIII along with the number of events sur- 
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TABLE 11: Numbers of observed events in data and expected 
yields from SM background processes for the twelve St and 
Ht bins in the ejj. channel. The quoted uncertainties are 
statistical only. 



St (GeV) 







0-70 70- 


-120 


>120 


(GeV) 


Data 


SM Data 


SM 


Data SM 


0-15 


1 


0.3 ±0.3 15 


13 ±2 


12 19 ±2 


15-60 


1 


0.09 ±0.1 6 


4.2 ±0.9 


11 8± 1 


60-120 





0.06±0.1 1 


1.6 ±0.6 


8 9±1 


>120 





0.01 ±0.05 


0.9 ±0.4 


6 7±1 



background, between 2% and 5% for the signal). These 
systematic uncertainties (except those for the luminos- 
ity and the instrumental background) are obtained by 
varying sequentially, before any selection, each concerned 
quantity within one standard deviation. For each chan- 
nel, the systematic uncertainty on the instrumental back- 
ground is estimated by varying the fit parameters within 
one standard deviation of their uncertainty. Higher sys- 
tematic uncertainties are observed for signal samples with 
low mi_^ and low Ato which give rise to higher event pop- 
ulations at low values of the px of the leptons and h 
quarks. 



viving at each step for the data, for each background 
component, and for the two signal samples A and B. 
Compared to the e/z channel, the estimated yields of 
tt, Z /^* —f T^T~ and diboson backgrounds are lower at 
the preselection stage. This is explained mainly by the 
threshold values of pr and ry used to identify electrons 
and muons. A slight excess of observed events is seen at 
the preselection level and is due to Z/^* e+e" events 
having no jets and for which the boson transverse mo- 
mentum is lower than 20 GeV. For these events, the par- 
ton showering implemented in the MC generators used in 
this analysis gives inaccurate results. The ti background 
dominates in the final stage of the selection. Four bins in 
Ht and St [Figs, [^e) and [IJf ) and Table |TV] are con- 
sidered to separate the signal from the SM background. 

For both e/i and ee channels, signal efficiencies, defined 
with respect to the numbers of events in the relevant 
channels, reach a value of 10% for large mass differences 
but decrease to values lower than 0.1% for Am < 20 GeV. 

The expected numbers of background and signal events 
depend on several measurements and parametrizations 
which each introduce a systematic uncertainty. The 
main sources of uncertainty that are common to efi and 
ee channels and affect both the backgrounds and the sig- 
nal consist of: electron identification and reconstruction 
efficiency (5% for the background, between 2% and 10% 
for the signal), jet energy cahbration (3% for the back- 
ground, between 2% and 11% for the signal), jet identifi- 
cation efficiency and energy resolution (2% for the back- 
ground, between 3% and 17% for the signal), luminosity 
(6.1%) [3], trigger efficiency (2%). The following sys- 
tematic uncertainties related to the background only are 
considered: instrumental background modeling (5% in 
the channel and 18% in the ee channel) and PDF 
(5% for diboson and 15% for tt and Z/j* processes). 
In addition, the e/i channel is affected by a systematic 
uncertainty related to the muon identification and recon- 
struction efficiency (2% for the background, between 2% 
and 5% for the signal). In the ee channel, an uncer- 
tainty coming from HF tagging is applied (2% for the 




Stop mass (GeV) 



FIG. 3: The 95% C.L. exclusion contour in the sneutrino mass 
versus stop mass plane. Shaded areas represent the kinemat- 
ically forbidden region and the LEP I [23j and LEP II Q ex- 
clusions. The dashed and continuous lines represent, respec- 
tively, the expected and observed 95% C.L. exclusion limit 
for this analysis. The band surrounding the observed limit 
denotes the effect of the uncertainty on the stop production 
cross section. 

No evidence for ii production is observed after ap- 
plying all selections for the e/x and ee data sets. No 
overlap is expected or observed between the two sam- 
ples. We combine the numbers of expected signal and 
background events and their corresponding uncertain- 
ties, and the number of observed events in data from 
the twelve bins of the e/i channel (Table [IT)) and the four 
bins of the ee channel (Table ITV|) to calculate upper lim- 
its on the signal production cross section at the 95% C.L. 
for various signal points using the modified frequentist 
approach This method employs a likelihood-ratio 
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TABLE III: Numbers of events observed in data and expected from SM background processes and the two signal samples A 
and B at the various stages of the analysis in the ee channel. The quoted uncertainties are statistical only. 



Selection 


Data 


Total SM 
Background Zj^* e+e~ 


Background contributions 
Z/7* T^T~ tt Diboson Instrumental 


Signal A 


Signal B 


Preselection 


27757 


25419 ±87 


24810 


120 


14.1 


23.4 


452 


10.7±0.5 


12.7±0.3 


Dielec 1 


6278 


6335 ± 38 


6143 


29 


14.2 


12.6 


136 


4.8 ±0.4 


10.6 ±0.3 


Dielec 2 


192 


200 ±5 


166 


11 


12.1 


3.9 


12 


3.0 ±0.3 


8.9 ±0.2 


Dielec 3 


142 


152 ±4 


122 


9.3 


11.4 


3.5 


5.8 


2.6 ±0.3 


8.0 ±0.2 


Dielec 4 


15 


16.0 ±0.6 


6.7 


0.5 


8.4 


0.22 


0.17 


0.6±0.1 


4.7±0.2 


Dielec 5 


12 


12.2 ±0.4 


3.0 


0.5 


8.4 


0.12 


0.16 


0.6±0.1 


4.6 ±0.2 



TABLE IV: Numbers of observed events in data and expected 
yields from SM background processes for the four St and Ht 
bins in the ee channel. The quoted uncertainties are statistical 
only. 



■5'T(GeV) 

Ht 45-150 >150 

(GeV) Data SM Data SM 
15-60 1 1.9 ±0.3 2 1±0.1 
>60 3 3.3 ±0.2 6 6 ±0.2 



(LLR) test-statistic, computed under the background- 
only (LLRb) or signal plus background (LLRg+b) hy- 
potheses. Simulated pseudo-experiments assuming Pois- 
son statistics and including the effect of systematic un- 
certainties are generated and distributions for LLRb and 
LLRs+b are obtained. By integrating the corresponding 
LLR distributions up to the LLR value observed in data, 
confidence levels CLb and CLg+b are derived. The stop 
cross section is varied until the ratio CL^, = CLg+b/CLb 
equals 0.05, which defines the 95% C.L. upper limit for 
the cross section for a given [m^^, mp] point. The in- 
tersection of the obtained cross section limit with the 
theoretical prediction for the cross section as a function 
of and yields the corresponding exclusion point 
in the [m^^ , nic:] plane. In this calculation, all systematic 
uncertainties except the ones related to the instrumen- 
tal background modeling and the PDF are considered 
as fully correlated between signal and background. The 
theoretical uncertainty of the stop signal cross section 
AfT^^I^ is estimated by adding in quadrature the vari- 
ations corresponding to the PDF uncertainty and the 
change in renormalization and factorization scale by a 
factor of two around the nominal value. Limits are es- 
timated for nominal (ar r ), minimal (err ~ -Aar ~ ) and 

maximal (ar i +AtTr r ) cross section values. We choose 
^ f iti . . 

not to correlate uncertainties between signal and back- 
ground so that the cross section limits can also be applied 
to other models or calculations. 

Figure [3] shows the excluded region as a function of 
the scalar top quark and sneutrino masses, for nomi- 
nal (continuous line) and for both minimal and maximal 



(band surrounding the line) values of a^^j^ , correspond- 
ing to the estimated theoretical uncertainty. For larger 
mass differences between the stop and the sneutrino, a 
stop mass lower than 175 GeV is excluded. A sensitiv- 
ity up to Am = 60 GeV is observed for stop masses of 
150 GeV. Combining the search in the ee final state with 
the e^ channel extends the final sensitivity by approxi- 
mately 5 GeV for large mass differences. The observed 
limit is within one standard deviation of the expected 
limit for m^^ > 150 GeV and within two standard devi- 
ations for TOj^ < 150 GeV. 

In summary, we presented the results of a search for 
the pair production of the lightest scalar top quark which 
decays into Mi'. Events with an electron and a muon or 
with two electrons have been considered in this analysis. 
No evidence for the lightest stop is observed in this decay, 
leading to a 95% C.L. exclusion in the [mi^,mc,] plane. 
The largest stop mass excluded is 175 GeV for a sneutrino 
mass of 45 GeV, and the largest sneutrino mass excluded 
is 96 GeV for a stop mass of 140 GeV. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 
(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina); FOM (The Netherlands); STFC (United 
Kingdom); MSMT and GAGR (Czech Republic); CRC 
Program, CFI, NSERC and WestGrid Project (Canada); 
BMBF and DFG (Germany); SFI (Ireland); The Swedish 
Research Council (Sweden); CAS and CNSF (China); 
and the Alexander von Humboldt Foundation (Ger- 
many). 



[a] Visitor from Augustana College, Sioux Falls, SD, USA. 

[b] Visitor from The University of Liverpool, Liverpool, UK. 

[c] Visitor from Rutgers University, Piscataway, NJ, USA. 

[d] Visitor from II. Physikalisches Institut, Georg- August- 
University, Gottingen, Germany. 

[e] Visitor from Centro de Investigacion en Computacion - 
IPN, Mexico City, Mexico. 



10 



[f] 
[g] 

[h] 
[i] 
It] 

[1] 
[2] 
[3] 
[4] 

[5] 

[6] 
[7] 
[8] 
[9] 

[10] 
[11] 



Visitor from ECFM, Universidad Autonoma de Sinaloa, 
Culiacan, Mexico. 

Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

Visitor from Universitat Bern, Bern, Switzerland. 
Visitor from Universitat Ziirich, Ziirich, Switzerland. 
Deceased. 

See, for instance, P. Fayet and S. Ferrara, Phys. Rep. 32, 
249 (1977); H.P. Nilles, Phys. Rep. 110, 1 (1984). 
See, for instance, A. Djouadi and Y. Mambrini, Phys. 
Rev. Lett. 63, 115005 (2001) and references therein. 
V.M. Abazov et al. (DO Collaboration), Phys. Lett. B 
659, 500 (2008). 

LEPSUSYWG Collaboration, ALEPH Collabora- 
tion, DELPHI Collaboration, L3 Collaboration, 
OPAL Collaboration, Note LEPSUSYWG/04-02.1, 



[12] 

[13] 
[14] 
[15] 
[16] 

[17] 
-418] 



'httpT/Tlepsusy. web . cern . ch/lepsusy / www/squarkslsumn^^ 
(DO Collaboration), Nucl. Instrum 



R. Hamberg, W.L. van Neerven, and T. Matsuura, 
Nucl. Phys. B359, 343 (1991) [Erratum- ibid. B364, 403 

(2002) ]. 

J.M. Campbell and R.K. Ellis, Phys. Rev. D 60, 113006 
(1999). 

N. Kidonakis and R. Vogt, Phys. Rev. D 68, 114014 

(2003) . 

R. Brun and F. Carminati, CERN Program Library Long 
Writeup W5013, 1993 (unpublished). 
The pseudorapidity rj is defined as r; = — ln[tan(S/2)], 
with 9 being the polar angle with respect to the 
proton beam direction. The rapidity y is defined as 
y = i ln[{E + pz)/iE-pz)]. 

G.C. Blazey et ah, in Proceedings of the Workshop: QCD 
and Weak Boson Physics in Run II, edited by U. Baur, 
R. K. Ellis, D. Zeppenfeld, Fermilab-Pub-00/297 (2000). 
V M Ahn^ov Pt al (DO Collaboration), Phys. Rev. D 74, 



V.M. Abazov et al. 
Methods Phys. Res. A 565, 463 (2006). 
W. Beenakker et al, Nucl. Phys. B515, 453 (1998). 
J. Pumplin et ah, JHEP 0207, 012 (2002). 
D. Stump et ah, JHEP 0310, 046 (2003). 
A. Pukhov et al, User's manual for version 3.3, INP-MSU 
98-41/542. 

T. Sjostrand et al, JHEP 05, 026 (2006); pythia version 
v6.319 is used. 

M.L. Mangano et al, JHEP 0307, 001 (2003); ALPGEN 
versions v2.05 and v2.11 are used. 



[19] T. Andeen et al, FERMILAB-TM-2365 (2007). 
[20] J. Smith, W.L. van Neerven, and J. A.M. Vermaseren, 

Phys. Rev. Lett. 50, 1738 (1983). 
[21] T. Scanlon Ph.D. thesis. Imperial College London, 

FERMILAB-THESIS-2006-43 (2006). 
[22] T. Junk, Nucl. Instrum. Methods Phys. Res. A 434, 435 
(1999); A. Read, in "1st Workshop on Confidence Lim- 
its," CERN Report No. CERN-2000-005 (2000). 
[23] C. Amsler et al, Phys. Lett. B 667, 1 (2008). 



